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ABSTRACT 

We use current measurements of the expansion rate H(z) and cosmic background 
radiation bounds on the spatial curvature of the Universe to impose cosmological 
model-independent constraints on cosmic opacity. To perform our analyses, we com- 
pare opacity-free distance modulus from H(z) data with those from two supernovae la 
compilations, namely, the Union2 and Sloan Digital Sky Survey samples. The influence 
of different SNe la light-curve fitters (SALT2 and MLCS2K2) on the results is also 
discussed. We find that these fitters present a significant conflict, with the MLCS2K2 
method being completely incompatible with a flat and transparent universe. 

Key words: cosmology - Distance scale - duality relation - galaxy clusters. 



> 

Ov 

t> 
o 

<N 



X 



1 INTRODUCTION 

Type la supernovae (SNe la) observations provide the most 
direct evidence for the current cosmic acceleration. In the 
context of Einstein's general theory of relativity, this result 
implies either the existence of some sort of dark energy, con- 
stant or varying slowly with time and space (see Caldwell 
and Kamionkowski (2010) and Li et al 2011 for recent re- 
views), or that the matter content of the universe is subject 
to dissipative processes (Lima & Alcaniz, 1999; Chimento et 
al 2003). 

However, there are still some possible loopholes in cur- 
rent SNe la observations and alternatives mechanisms con- 
tributing to the acceleration evidence or even mimicking the 
dark energy behaviour have been proposed. Examples are 
possible evolutionary effects in SNe la events (Drell, Loredo 
& Wasserman 2000; Combes 2004), local Hubble bubble 
(Zehavi et al. 1998; Conley et al. 2007), modified gravity 
(Ishak, Upadhye k Spergel 2006; Kunz k Sapone 2007; 
Bertschinger & Zukin 2008), unclustetered sources of light 
attenuation (Aguirre 1999; Rowan- Robinson 2002; Goobar, 
Bergstrom & Mortsell 2002) and the existence of Axion-Like- 
Particles (ALPs), arising in a wide range of well- motivated 
high-energy physics scenarios, and that could lead to the 
dimming of SNe la brightness (Avgoustidis et al. 2010). 

On the other hand, several authors have recently dis- 
cussed how the so-called cosmic distance duality (CDD) 
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relating the luminosity distance (Dl) to the angular diam- 
eter distance (Da) of a given source can be used to verify 
the existence of exotic physics as well as the presence of sys- 
tematic errors in SNe la observations (Basset & kunz 2004; 
Uzan, Aghanim & Mellier 2004; Holanda, Goncalves & Al- 
caniz 2012). As is well known, the CDD relation is closely 
connected with the Etherington's reciprocity law (Ethering- 
ton 1933, Ellis 2007), being valid for all cosmological mod- 
els based on Riemannian geometry and requiring only that 
source and observer are connected by null geodesies in a 
Riemannian spacetime and that the number of photons is 
conserved (Ellis 2007). 

The implementation of the CDD tests mentioned above 
follows different routes. Basset & Kunz (2004), for instance, 
used SNe la data as measurements of Dl and estimates of 
Da from FRIIb radio galaxies (Daly & Djorgovski 2003) and 
ultra compact radio sources (Gurvitz 1994; Lima & Alcaniz 
2002; Alcaniz 2002) to test possible new physics with ba- 
sis on a generalization of the CDD relation. Uzan, Aghanim 
& Mellier (2004) argued that the Sunyaev-Zel'dovich effect 
plus X-ray techniques for measuring Da from galaxy clus- 
ters is strongly dependent on the validity of this relation. By 
assuming a deformed CDD relation, Dl(1 + z)~ 2 /Da = rj, 
in a ACDM background they found that the value rj = 1 
is only marginally consistent with the galaxy clusters data 
used in the analysis (for the role of the cluster geometry in 
this test, see Holanda, Lima & Ribeiro, 2010; 2012). More 
recently, a cosmological model-independent test involving 
only measurements of the gas mass fraction of galaxy clus- 
ters from Sunyaev-Zeldovich and X-ray surface brightness 
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observations was discussed by Holanda, Goncalves & Al- 
caniz (2012). From this analysis, no significant violation of 
the CDD relation was found, with the value 77 = 1 lying 
in the la interval (for other CDD analyses, we refer the 
reader to Holanda, Lima & Ribeiro, 2011; Nair, Jhingan & 
Jain 2011; Li, Wu & Yu 2011; Meng et al. 2011; Goncalves, 
Holanda & Alcaniz, 2011; Lima, Cunha & Zanchin, 2011). 

In a different approach, Avgoustidis et al. (2010) ex- 
plored consistency among different distance measurements 
by considering a possible violation of cosmological photon 
conservation as the only source of CDD violation. The au- 
thors assume a flat ACDM model and the distance parame- 
terization, Dl — Da(1 + ^) 2+c , to place bounds on the CDD 
parameter e by combining the SNe la Union compilation 
(Kowalski et al. 2008) with the latest measurements of the 
Hubble expansion in the range < z < 2 (Stern et al. 2010). 
The basic idea behind the test proposed is that, while SNe 
la observations are affected by at least four different sources 
of opacity, namely, the Milky Way, the hosting galaxy, in- 
tervening galaxies, and the Intergalactic Medium, the cur- 
rent H(z) measurements are obtained from ages estimates 
of old passively evolving galaxies, which relies only on the 
detailed shape of the galaxy spectra, not on the galaxy lu- 
minosity. Therefore, differently from Dl measurements from 
SNe la, H(z) observations are not affected by cosmic opac- 
ity t(z) since this quantity is assumed to be not strongly 
wavelength dependent on the optical band (see Avgoustidis 
et al. 2009 and references therein for more details). By fit- 
ting the ACDM model to H(z) data and using the best-fit 
values to estimate Dl, a direct comparison is made with 
SNe la observations to constrain the cosmic opacity. From 
a direct relation between the CDD parameter e and r(z), 
Avgoustidis et al. (2010) found e = -OM±° ° 8 7 (2a). 

In this paper, we discuss a cosmological model- 
independent version of the CDD test proposed by Avgous- 
tidis et al. (2010). Our approach differs from the one de- 
scribed above in that opacity- free estimates of Dl are ob- 
tained from a numerical integration of current H(z) data 
points and not from a fit of a given cosmological model to 
the data. Constraints on e are derived by comparing these 
estimates with SNe la observations from the Union2 sample 
(Amanullah et al. 2008) and two compilations of the nearby 
+ SDSS + ESSENCE + SNLS + Hubble Space Telescope 
set (Kessler et al. 2009) - with which we discuss the influ- 
ence of the different SNe la light-curve fitters (SALT2 and 
MLCS2K2) in the results. From the analyses performed, we 
find that evidence for cosmic opacity depend strongly upon 
the SNe la light-curve fitting, with SNe la samples that use 
SALT2 prefering a transparent universe and the MLCS2K2 
compilation ruling out such a possibility by ~ 3.5cr. 



LUMINOSITY DISTANCE FROM H(Z) 
MEASUREMENTS 




Figure 1. Dl(z) obtained from measurements of the expansion 
rate (filled circle). For the sake of comparison, we also plot Dl 
measurements extracted from the Union2 SNe la sample (open 
circles). The curves stand for the second degree polynomial fit of 
the opacity-free Dl points and the corresponding lcr error. 



archival data (Abraham et al. 2004; McCarthy et al. 2004), 
Stern et al. (2010) have calculated 12 H(z) data points over 
the redshift range 0.1 < z < 1.75 0- 

In order to transform these H(z) measurements into dis- 
tance estimates, we solve numerically the comoving distance 
integral for non-uniformly spaced data 
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using a simple trapezoidal rule. Since the error on z mea- 
surements is negligible, we only take into account the un- 
certainty on the values of H(z). As one may check, by using 
standard error propagation techniques, the error associated 
to the i th bin is given by 
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so that the error of the integral J2} in the interval z = 
- z n is S n = Yl™-! Si. In the analyses perfomed, we added 
to our H(z) sample, the recent value of the expansion rate 
today Ho = 73.8 ± 2.4 km/s/Mpc, as obtained by Riess et 
al. (2011). The Dl data points obtained from H(z) obser- 
vations are shown in Fig. 1. For the sake of comparison, we 
also plot Dl measurements extracted from the Union2 SNe 
la sample. 



Current measurements of the expansion rate at z 7^ are 
obtained by calculating the derivative of cosmic time with 
respect to redshift, i.e., H(z) ~ — (itp^y • This method was 
first presented by Jimenez and Loeb (2002) and consists in 
measuring the age difference between two red galaxies at 
different redshifts in order to obtain the rate Az/At. By 
using a recent released sample from Gemini Deep Survey and 



1 Direct measurements of H(z) at different redshifts will also be 
possible through measurements of the line-of-sight or radial com- 
ponent of baryonic acoustic oscillations (BAO) from large redshift 
surveys with redshift precision of the order of 0.003(1 + z) - see, 
e.g., Benitez et al. 2009. 
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Figure 2. Confidence contours on the plane Qj~ — e for the three SNe la sub-smaples discussed in the text. Contours are drawn for 
X 2 = 2.30 (2a) and 6.17 (2a). The point where the two solid lines crosses stand for a perfect transparent (e = 0) and flat universe 
(Qfc = 0). From these panels, it is clear that evidence for cosmic opacity depend strongly upon the SNe la light-curve fitting. 



3 CONSTRAINTS ON COSMIC OPACITY 
3.1 Methodology 

As argummented by Avgoustidis et al. (2010), if there were 
a source of photon absorption affecting the universe trans- 
parency, the distance modulus derived from Supernovae 
would be systematically affected. In particular, any effect 
that reduces the number of photons would dim the SNe la 
brightness and increases Dl- Thus, if r(z) denotes the opac- 
ity between an observer at z = and a source at z due to, 
e.g., extinction, the flux received from the source would be 
attenuated by a factor e~ T ^ and thus the observed lumi- 
nosity distance (DL,obs) is related with the true luminosity 
distance (D L ,true) by 

D 2 L,obs — D 2 L ^ true e T ^ . (4) 

Therefore, the observed distance modulus is given by (Chen 
and Kantowski, 2009a; 2009b) 

rriobs(z) = mtrue(z) + 2.5 [log e]r(z) . (5) 

In our analyses, measurements of m bs are taken from the 
SNe la Union2 and SDSS compilations (see Section 3.2). 
Note that, differently from the Avgoustidis et al. (2010) 
analysis, which fitted a flat ACDM model to H(z) data 
and used the best-fit values to build a DL,tme sample, we 
compare DL,obs estimates from SNe la data to opacity- free 
luminosity distance DL,true inferred directly from the H(z) 
measurements, as described earlier. Note also that, since the 
SNe la and H(z) observations are performed at different z, 
we calculate DL,true at each SNe la redshift from a second 
degree polynomial fit of the data points shown in Fig. 1. 

In order to proceed further, we must assume an appro- 
priated redshift dependence for r(z). Here, we follow Av- 
goustidis et al. (2009; 2010) and consider 

r(z) = 2ez . (6) 

As shown in the above reference, this linear expression can 
be derived from the usual CDD parameterization Dl = 
Da(^ + z) 2+e for small values of e and z < 1, where e quanti- 
fies departures from transparency. To discuss our results in a 



more general setting, we also consider Fridmann- Robertson- 
Walker geometries with arbitrary curvature and use in our 
analyses the current WMAP bound on the curvature param- 
eter, i.e., Q k = 0.0125 ± 0.0155 (2a). 



3.2 Data sets and analysis 

As mentioned earlier, we use two SNe la data sets in our 
analyses. The Union2 sample is an update of the original 
Union compilation that comprises 557 data points including 
recent large samples from other surveys and uses SALT2 for 
SNe la light-curve fitting. The second sample is the nearby 
+ SDSS + ESSENCE + SNLS + Hubble Space Telescope 
set of 288 SNe la (throughout this paper we refer to this 
set as SDSS compilation) that uses both SALT2 (Guy et al. 
2007) and MLCS2K2 ( Jha et al. 2007) light-curve fitters 
and is distributed in redshift interval 0.02 ^ z < 1.55. In 
order to verify the effect of the light-curve fitting on the 
results, we consider both SDSS sub-samples in our analyses. 

We estimated the best-fit to the set of parameters P = 
(Qfc,e), by evaluating the likelihood distribution function, 
C oc e~ x / 2 , with 

2 _ (m obs (z) - mtrue(z) - 2.17ez) 2 

X ~ 2-^ (j 2 + (T 2 + (T 2 ' ^ 

z U m(obs) ^ U m(true) ^ U int 

where cr^ (true) , cr^ (obs) and a 2 nt are the errors associated 
to distance modulus from H(z) measurements, distance 
modulus from SNe la and additional intrinsic dispersion 
to each SNe la sample, respectively. m tr ue is obtained via 
TTltrue — 51og 10 D L ,true + 25, while D L ,true is given by one 
the following forms: 



2 MLCS2K2 calibration uses a nearby training set of SNe la as- 
suming a close to linear Hubble law, while SALT2 uses the whole 
data set to calibrate empirical light curve parameters, and a cos- 
mological model must be assumed in this method. Typically a 
ACDM or a c^CDM model is assumed. Consequently, the SNe la 
distance moduli obtained with SALT2 fitter retain a degree of 
model dependence (Bengochea 2011). 



© 2011 RAS, MNRAS 000, [T]-?? 



4 R. F. L. Holanda, J. C. Carvalho and J. S. Alcaniz 




Figure 3. Likelihood function for e after marginalizing over the curvature parameter. 
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where = cH 1 and Dc was denned in Eq. J2]). 



3.3 Results 

The results of our statistical analyses are shown in Figures 2 
and 3. Figures 2a-2c show contours of A% 2 = 2.30 and 6.17 
on the Qk - e plane when the Union2, SDSS (SALT2) and 
SDSS (MLCS2K2) compilations are considered, respectively. 
For the Union2 sample, we find that a perfect transparent 
(e = 0) and flat universe is allowed by the current data 
at 2a level, with e = 0.03 ± 0.02 (la). The influence of 
the light-curve fitting on the analysis can be easily seen by 
comparing Panels 2b and 2c. Note that, while the SDSS 
sample that uses SALT2 is in full agreement with a flat 
universe with perfect cosmic transparency, the MLCS2K2 
compilation rules out such a possibility by ~ 3.5cr. For these 
two analyses, we find e = 0.025 ±0.023 and e = 0.105 ±0.03 
at 68.3% C.L., respectively, which clearly shows conflicting 
results between the two light-curve fitters. 

By marginalizing over the curvature parameter, Pan- 
els 3a-3b display the likelihood for the e parameter for 
Union2, SDSS (SALT2) and SDSS (MLCS2K2) compila- 
tions, respectively. In this case we obtain e = 0.030 ±0.015, 
e = 0.025±0.021 and e = 0.105±0.025 at la level. It is worth 
mentioning that these bounds on e are only marginally com- 
patible with those found by Avgoustidis et al. (2010), i.e., 
e = — 0.040 Q7 at 2a level. For the sake of comparison, we 
also repeat our analysis using the Union compilation (the 
same used by Avgoustidis et al. 2010). For this analysis, we 
obtain e = 0.013 ± 0.013 at la level. 



4 CONCLUSIONS 

The results of observational cosmology in the last years have 
opened up an unprecedented opportunity to test the verac- 
ity of a number of cosmological theories as well as the ex- 
istence of new physics in the Universe. Motivated by these 



results, several analyses have recently discussed the impor- 
tance of comparing cosmological distances to explore devi- 
ations from the standard cosmological model or a possible 
presence of exotic physics. However, most of these analyses 
are cosmological model-dependent, which makes the results 
less general than desired. 

In this paper, we have used recent H(z) measure- 
ments from passively evolving galaxies to obtain cosmolog- 
ical model-independent distance modulus and impose con- 
straints on cosmic opacity by comparing these data with 
current SNe la observations. In order to perform our anal- 
ysis, we have considered two of the most recent samples of 
SNe la, namely, the Union2 and the SDSS compilations. We 
have found that both Union2 and SDSS (SALT2) compila- 
tions are in full agreement with a perfect transparent and flat 
universe whereas the SDSS compilation that uses MLCS2K2 
light-curve fitting rules out such a possibility by ~ 3.5cr. If 
upcoming SNe la and H(z) data confirms the MLCS2K2 
results, this would open room for existence of exotic physics 
as suggested by SNe la observations. In this regard, a pos- 
sible coupling of photons to particles beyond the standard 
model of particle physics modifying the apparent luminosity 
of sources has recently been considered (for a recent review 
on this subject see, e.g., Jaeckel and Ringwald, 2010). 
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